Jet rates in the hard scattering process at finite temperature 

Krishnendu MukherjeeQ 
Department of theoretical Physics, Tata Institute of Fundamental Research, Homi Bhabha Road, 

--^- . Mumbai 400 005, India. 

O 

o 

(N 
Q^' Abstract 

(D 

We compute the cross-section of the hadronic jets arising from the quark antiquark pair which 
are produced from a hard photons (of 4-momentum q) in the plasma, predominantly consisting of 
K»" ! thermalised quarks and gluons. The quark antiquark pair is hard and scattered off the heat bath to 

^ \ form jets, while the gluons being soft get thermalised in the heat bath. The infrared divergences cancel 

Q ■ in the observable cross-section to as order, which includes the process of emission and absorption of 

Q . real gluons. Since the massless quark antiquark pair is hard the Compton scattering processes are 

i-^ \ absent in the heat bath and it renders an uncancelled collinear divergent piece in the cross-section. 

Q_i' We regularize it by eliminating the collinear region from the phase space and write it in terms of 

^ ■ jet parameters. The temperature dependent part of the jet cross-section is regular at large -^ and 



vanishes when -^ > oo. Since jets carry the thermal signature of the hot plasma the jet production 

XI ■ 

;h ' rate can be used as a thermometer of the heat bath. 



In this paper we consider a heat bath containing plasma of quarks, antiquarks and gluons at 

temperature T (= 1//3). In the e^e~ collision a hard photon (hard means it's energy and momentum 

are much greater than T.) is produced. Photon enters into the heat bath and produce a hard quark 

antiquark pair and a soft gluon (soft means it's energy and momentum are much less than T.). Since 

the gluon is soft it gets thermalised, however the quark antiquark pair being hard, leave the heat bath 

without loosing much of their initial energy and form jets. We have used the real-time formalism 

of finite temperature field theory to calculate the imaginary part of the virtual photon polarization. 

We take quarks as massless. The cross-section for the real processes (7* — > qqg and 7*5' -^ qq) are 

infrared divergent. The divergences neatly cancel from the virtual processes (7* -^ qq) at a^ orderQ. 
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^At temperature T = the KLN Q] theorem ensures that the physically measurable probabilities are infrared finite 



However there still remains an uncancelled part, which is a collinear divergent piece, appearing from 
the real emission process (7* — > qqg). Since the quark antiquark pair is hard the Compton scattering 
processes (7*^ -^ qg and 7*g -^ qg) are absent in the plasma. As a result we are left with this 
collinear divergent piece in the real emission process and hence also in the cross-section. We have 
regularised this divergent piece by eliminating the collinear region from the phase space and write it 
in terms of the jet parameters. Since jets carry the information about the plasma in the heat bath, 
the jet production rate (or cross-section) can be used as a thermometer of the heat bath0 ||1 1|| . 

We consider the hard scattering of quark antiquark pair in the plasma. Since the gluons are very 
soft they get thermalised in the heat bath. In general, if both quarks and gluons are thermal the 
imaginary part of the photon self energy can be written as 0, ^ 

ImU,^{q) = 2^^w^ni^,(g) (1) 

where Uj^^, is the 12 component of the self energy. In order to obtain the imaginary part for the 
non-thermal quarks and the thermal gluons we will set all the sign function in the quark propagators 
to theta function and also set all the fermion distribution functions to zeroQ. 
The cross-section of e~^e~ — ;> 7* — ;> qqg is given as 

a = 2ImI[^,{q)L^\ (2) 

L'^'^ is the leptonic part and is given as 



L"^ = TTT^ (P1P2 + P1P2 - V,uPi ■P2), (3) 



;2 

2{q' 

where pi and p2 are the 4-momenta of e^ and e~ and q {= Pi + P2) is the 4-momentum of the 
virtual photon. The imaginary part of the photon self energy, Imll^^{q), contains the hadronic 
even when masses go to zero. Infrared divergences become severe at T 7^ than at T = because the phase space 
integrals are weighted by the Bose-Einstein (BE) distribution function. In order to get infrared cancellations at T 7^ 
KLN prescription is to include in the observable rate the emission and absorption of low energy real particles in the 

plasma [2-12]. 

^Temperature of the plasma can also be measured by studying the missing-p^ spectrum of multijet events where 

Pt is the transverse component of the total momentum of the jet ||13| . 

^ Massless fermion and gauge field (in Feynman gauge) propagators at finite temperature read iSAB{p) ~ i> iSAB{p) 

and iD^^ (p) — —■q'^'-'^ 5°"^ "iDab {A,B = 1,2) respectively, where t?'"' is the Minkowski metric and a,b {a,b = 

1, • • • , N^ — 1 for SU{Nc) gauge group) are the gauge group indices. iSab and iDAB are the 2x2 matrices and their 

explicit form can be found in Ref . |]l5|, ll^ . 



contributions. Since the electromagnetic current is conserved we can write it as 

Jmn,,(g) = -^ (^-r^,, + '^^ ImW^{q) (4) 

where c? (= 4 — 2e, e > 0) is the space-time dimension. Since the leptons are taken to be massless, 
the cross-section reads 

In the following we calculate the scalar ImU'^ to obtain the jet cross-section to ag order. We shall 
work here in the rest frame {u = 0) of the heat bath. 

The loop diagram corresponding to the Born amplitude of the process 7* -^ qq is shown in the 
Figl. The imaginary part of the trace of photon polarisation tensor read 



/mn(^j^(g) = {2-d)e'N,{J2Ql)q'ldR2 



^^Vf^V'IilziL^e (6) 



47r ^ V^TT/iV r(2-2e) 
where eQ„ is the charge of the n-th flavour quark and indices n runs over the number of quark 
flavours. dR2 is the differential of two particle phase space for massless quarks. Since the quarks are 
hard they are non-thermal and the trace of the polarisation tensor is the same as zero temperature 
result. 

The two-loop diagrams corresponding to the real processes 7* -^ qqg and '~f*g -^ qq are shown 
in the Fig2. The complete contribution of these real diagrams to the imaginary part of the trace of 
photon polarisation tensor is given as 

2 2 
Jmn(«^'^(g) = Y. ImU^^'l^iq) = ^J-C{R)N,{Y^ Ql) [j(°) + j(^™) + J^'^^^)] , (7) 

j=a,b,c,d n 

where 

J(°) = J dRsX , J^^'") = J dRsnB{Ks)X , J^^"'^ = J dRsnB{Ks)X, (8) 

and 

X{K,,K2,Ks) = -—^^—^-—\q\K^.K2) + {K^.K.f + {K^.K,f - e{K^.K, + K^.K.f]. (9) 

Here C{R) is the quadratic Casimir invariant of the representation R of the quarks. For SU{Nc) 
gauge group C{R) = {N^ — l)/2Nc, where the quark furnishes the fundamental representations of 
the group. Ki and K2 are the 4-momenta of the quark and antiquark respectively and K3 is the 



4-moinentum of the gluon. dR^ and dR^ are the differential 3-body phase space for emission and 
absorption processes respectively for massless quarks and gluonsQ. 
After performing the integrations the final results are 
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Virtual diagrams are shown in the Fig.3 and Fig.4. First we consider the vertex correction 

diagrams shown in the Fig.3. The complete contribution of the vertex correction diagrams to the 

imaginary part of the trace of photon polarisation tensor is given as 

= ImI\.^^^^^^'{q) + ImI\.^^^^^^'{q). (13) 

Here the ImYv^^^^ is written as the sum of the two terms: Zero temperature part and the finite 
temperature part. The expression for them are the following: 
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(14) 
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^ 3-body phase space for emission process has the momentum conserving delta function 6^'^'{q — Ki — K2 — K3) in 
d dimension. To evaluate the phase space integrals for emission precess one defines the new variables Xi — 2q.Ki/q^^ 
j = 1, 2, 3 where xi + 0:2 + 2^3 = 2 ||l^, Q. For absorption process the 3-body phase space integral contains 6^'^'>{q — 
Ki — K2 + K-^). To evaluate the 3-body phase space integral for absorption process one goes over to the new variable 
Xi {i — 1, 2, 3) where Xi now satisfy the relation xi + X2 — x^ — 2. 



Next consider the self energy correction diagrams in Fig.4. The complete contribution of the 
diagrams in the Fig.4 to the imaginary part of the trace of photon polarisation tensor is given as 



/mn<*'''(,) 



^c(fi)iv.(i:«),' 



g2 \-^' r2(2-e) 2(1 -e) 



ATTfiy r2(2-2e) e 



X / dXgXg ^'nBiJq^X3/2). 



(16) 



To obtain the total cross-section for the processes e~^e~ — ;> 7* — ;> qqg and e~^e~ -^ '-y*g -^ qq we 
have to sum the real and the virtual process contributions. Therefore using the formula in eqn.(^) 
the total cross-section to a^ order reads 

^2 3 rl 
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where a = j-, as = §- and cr^ is called the Born cross-section and is given as 
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Here F contains the finite temperature divergent contributions from the real and the virtual 
diagrams. Its complete expression is 

■ ^2 \ - rfi - 
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After the cancellation of infrared divergences between the virtual and real processes finite part of F 
reads 
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All the terms in eqn.(^) are finite except for the first one. The integral in the first term is finite, 
however due to the - factor before it, it becomes divergent when e ^ 0. This divergence appears 
when the gluon momentum is parallel to quark or antiquarkQ. The natural way of regularizing this 
divergence is to eliminate the potentially dangerous collinear integral region from the phase space. 
One can eliminate the regions around 6 = {0 < 6 < 9.^) and 6 = -k {-k — 9^ < 9 < -k) in the 
angular integral. Then the infrared divergence would have manifested itself as singularity in 9m for 
its vanishing limit. If the resolving power of the detector is less than 9m it can not be able to detect 
the soft gluons as separate entity from the quarks. So we can relate 9m to the maximum power of 
the detector to detect the two particles separately, which are coming towards it with a small angular 
separation between themselves. In order to express this - factor in terms of 9m we look back into 
the evaluation of j('^'") and the relation is 
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Then the cross-section reads 

a = a(°)+a(^), (22) 

where a'^'^^ and a'^'^'^ are the zero temperature and finite temperature piece of the total cross-section 
respectively. The expression for them are 
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^This is a collinear divergent piece which appears due to the absence of Compton scattering processes in the heat 
bath 



We are considering here the hardonic jets arising from the quark-antiquark pair in e"'"e~ annihila- 
tion. In order that the hadronic jets follow from the quark-antiquark pair, the quark and antiquark 
are required not to loose too much energy in their direction by the emission of gluons and quark 
pairs. In otherwords most of the annihilation energy is deposited along the direction of the quark 
and antiquark. For the production of 2-jets the three body phase space integration in eqn.(|^) for 
emission process must be restricted. We define the contribution of the real diagrams to the imaginary 
part of the trace of photon polarisation tensor as 

/^ng;-(g) = '^CiR)N^{Y: Ql) [4^ + Ji"^^ + J("^^)] (25) 

^ n 

where 

4°) = / dR^X, (26) 

Jn 

j{em) 



jr' = / dR,nBiK,)X, (27) 

X is given by eqn.(y) and the integral region TZ is specified by the following conditions: The emitted 
gluon is either soft (i.e. 0:3 < A) or collinear to one of the quarks (i.e. 612, O23 < 26) where 0:3 
(= 2K3.q/q^) is the fraction of the energy of the virtual photon carried by the gluon and ^13 (^23) is 
the angle between the gluon and the quark (antiquark). Evaluating the integrals we obtainQ 
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^ In addition to Xi {— 2q.Ki/q^ where i = 1,2, 3) define two more variables ijlj 6 — (1 — cos 9iz)/2 (/ ~ 1,2). Among 
these five variables two are independent and so we choose ^1 and X3 as the independent variables. Region TZ is specified 
by imax < 6 < imm and < 0:3 < A, where £^,nax = i_^3^(^I'"3/sin:jj and S^^mn = sin^ 5. Now J^^ = j(°) - J^*', where 
iZ is obtained by eliminating TZ from the whole phase space. So J^ = iq{2ttY^ /a '^^'i^ii^ ^ 2:3)/,'"°^ d^i (i_^ ^ -.a 
where X is written in terms of 2:3 and ^1. Evaluation of J^j" is also carried out in the same fashion. 



The 2-jets cross-section for the processes e'^e 
reads 
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qqg and e~^e -^ 7*^7 — > qq to a^ order 
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where a2jit and o"2jei are the zero and finite temperature piece respectively. The expression for them 



are the following: 



cr. 



a^ 



(0) 
2jet 

(T) 
2jet 



O-B 



1 + — C(i?) 7-^^ + --41n(51nA-31n(5 
TT 16 2 



a.. 



as— C(i?) 

vr 



tvr^^ + /2(v/^/r,A,<5) 



(31) 
(32) 



where 
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and (Tb is given by the eqn. (p!8[) . Here without any loss of generality we have chosen 6m = 2(5. In 
real life we can relate 6 to the maximum resolving power of the detector. According to eqn. (^) and 
eqn. (|32D we see that the order a^ correction to the two jets from the quark pair is controllable in 
size within the frame work of perturbation theory at finite temperature. Since the virtual photon is 



hard, the ratio -^7^ is very large. The function f2iVq^/T) is well behaved for large -Kjt- ratio and it 



vanishes for 



00. 



3-jets cross-section is the complementary of the 2-jets cross-section. For 3-jets the imaginary part 
of the trace of photon polarisation tensor reads 
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The 3-jets cross-section for the process e'^e — > 7* — > qqg to as order reads 

= 4;l+^s, (35) 

where a^^jl^ and o"3jgi are the zero and finite temperature piece respectively. The expression for them 



are 
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where 



/3(^/T) = - /' ^ns{^'xs/2){xl + 2((1 - 0:3)^ + 1) ln5 + ((1 - x^f + 1) ln(l - x-M. (38) 



Here fz{\f(p /T) is also well behaved for large ^^^ ratio and it vanishes for ^ > 00. Because of 

In 5 and ln(l —0:3) in the second and third term within the second bracket of eqn. (|38|) respectively we 
obtain positive contribution after integrations. However the first term within the bracket will give 
negative contribution. We write the 3-jets cross-section as 

assets = aB—CiR)h{^qyT, 6, A) + 0(5, A), (39) 

where 

;i = 41n51nA + 31n5+^-^ + /3(V^/T). (40) 

The plots of h versus 6 (Fig. 5 and Fig. 6) for different set of values of \/q^/T and A show that the 
function h remains positive for realistic values of its argument. The maximum value of A we have 
taken is 0.4, which is much away from the soft gluon approximation, yet h remains positive for this 
large value of A. 

We obtain from eqn. (^) and eqn. (|35|) that the ratio of the two cross-section is 

^m = ^ciRM^yr, 5, A) + oial). (4i) 

0'2jet TT 

Apart from the small numerical factor ^C{R) this ratio behaves in the same way as the function 
h. As 6 and A decrease the 3-jets phase space gets increased compared to the 2-jets, as a result this 
ratio gets increased which are evident from the plots in Fig.5 and Fig. 6. It is also evident from these 



plots that for a fixed value of 6 and A as -y/g^/T increases the phase space of gluons get decreased 
due to Bose-Einstein distribution function and as a result this ratio gets decreased. 

I gratefully acknowledge S. Gupta for suggesting the problem and for the many subsequent 
discussions and helpful comments. I also thank P. Roy and K. Sridhar for helpful discussions and 
suggestions. I am also thankful to P. Aurenche for clarifying certain points regarding the infrared 
cancellations in their paper P] and also for helpful comments. 

References 

[1] T. Kinoshita, J. Math. Phys. 3, 650 (1962); T. D. Lee and M. N. Nauenberg, Phys. Rev. 133, 
B1549 (1964). 

[2] R. Baier, B. Pire and D. Schiff, Phys. Rev. D 38, 2814 (1988). 

[3] T. Altherr, P. Aurenche and T. Becherrawy, Nucl. Phys. B315, 436 (1989). 

[4] T. Altherr and P. Aurenche, Z. Phys. C 45, 99 (1990); T. Altherr and T. Becherrawy, Nucl. 
Phys. B330, 174 (1990); J. Cleymans and I. Dadic, Z. Phys. C 42, 133 (1989); Y. Gebellini, T. 
Grandou, D. Poizat, Ann. Phys. (N.Y.) 202, 436 (1990). 

[5] T. Grandou, M. Le Bellac and J. L. Meunier, Z. Phys. C 43, 575 (1989); T. Grandou, M. Le 
Bellac and D. Poizat, Phys. Lett. B 249, 478 (1990); Nucl. Phys. B358, 408 (1991). 

[6] E. Braaten and R. D. Pisarski, Phys. Rev. Lett. 64, 1338 (1990); Nucl. Phys. B337, 569 (1990); 
B399, 310 (1990). 

[7] E. Braaten and M. H. Thoma, Phys. Rev. D 44, 1298 (1991); 44, 2625 (1991). 

[8] T. Altherr, Phys. Lett. B 262, 314 (1991). 

[9] H. A. Weldon, Phys. Rev. D 44 3955 (1991). 
[10] P. V. Landshoff and M. Taylor, Nucl. Phys. B430, 683 (1994). 

[11] S. Gupta, D. Indumathi, P. Mathews and V. Ravindran, Nucl. Phys. B458, 189 (1996). 
[12] D. Indumathi, Ann. Phys. (N.Y.) 263, 310 (1998). 



[13] S. Gupta, Phys. Lett. B 347, 381 (1995). 

[14] R. L. Kobes and G. W. Semenoff, Nucl. Phys. B260, 714 (1985); Y. Fujimoto, M. Morikawa 
and M. Sasaki, Phys. Rev. D 33, 590 (1986). 

[15] For a review of finite temperature field theory, see N. P. Landsman and Ch. G. van Weert, Phys. 
Rep. 145, 141 (1987). 

[16] A. J. Niemi and G. W. Semenoff, Ann. Phys. (N.Y.) 152, 105 (1984); Nucl. Phys. B230, 189 

(1984). 

[17] T. Muta, Foundations of Quantum Chromodynamics (World Scientific Lectures Notes in Physics 
- Vol. 5, 1987). 

[18] R. D. Field, Applications of perturbative QCD (Addision- Wesley, 1989). 




Figure 1: One loop photon polarization graph for the computation of Born amplitude. 




(c) 




Figure 2: Two loop photon polarization graph for the computation of'j* — * qqg and '~f*g —>■ qq 

amplitude . 





Figure 3: Two loop photon polarization graph for the computation of vertex correction of'~f* -^ qq 

amplitude to a^ order. 






(b) 




Figure 4: Two loop photon polarization graph for the computation of self energy corrections of 

7* -^ qq amplitude to as order. 
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Figure 5: Plot of h versus 6 for different values of A and ^fifjT = 1.2. 
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Figure 6: Plot of h versus 6 for different values of A and \f(^jT = 10. 



